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ABSTRACT 



The recently discovered exoplanet G1581d is extremely close to the outer edge of its system's habitable zone, which has led to 
much speculation on its possible climate. We have performed a range of simulations to assess whether, given simple combinations 
of chemically stable greenhouse gases, the planet could sustain liquid water on its surface. For best estimates of the surface gravity, 
surface albedo and cloud coverage, we find that less than 10 bars of CO2 is sufficient to maintain a global mean temperature above the 
melting point of water. Furthermore, even with the most conservative choices of these parameters, we calculate temperatures above 
the water melting point for CO2 partial pressures greater than about 40 bar However, we note that as G1581d is probably in a tidally 
resonant orbit, further simulations in 3D are required to test whether such atmospheric conditions are stable against the collapse of 
CO2 on the surface. 
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1. Introduction 

In 2007, radial velocity measurements were used to discover two 
new planets in the G1581 system (Udry et al. 2007). These plan- 
ets have captured much attention both in the community and 
among the general public, as their minimum masses were mea- 
sured to be below 10 MEanh, and they are close to the edges of 
their system's nominal 'habitable zone', i.e., the loosely defined 
orbital region in which planets can sustain liquid water on their 
surfaces. The first planet, G1581c, which is closer to its star and 
was the first discovered, was initially estimated to be potentially 
habitable based on its equilibrium temperature T^q - 320 K, us- 
ing an Earth-like planetary albedo 0.29. In contrast, the second 
planet G1581d has an equilibrium temperature T^q = 195 K for 
an albedo of 0.2, which suggests it may be too cold to sustain 
surface liquid water. However, these analyses neglect any possi- 
ble warming of the surface due to the planet's atmosphere. 

In the first detailed assessment of the potential habitability 
of these planets, [Selsis et aL] ( |2007| l reviewed a variety of factors 
that could influence their climates. They concluded that based 
on standard assumptions of atmospheric warming by a mix- 
ture of CO2 and H2O (with possible regulation of CO2 via the 
carbonate-silicate cycle), G1581c was unlikely to be habitable, 
while for G1581d, the situation was much less clear. According 
to Kasting et al. ( 1993| ), the outer edge of the habitability zone 
is most likely the distance at which CO2 condensation begins to 
occur on the surface of the planet. However, CO2 condensation 
in the atmosphere leads to the formation of CO2 clouds, which 
can cause a strong warming effect due to the scattering of in- 
frared radiation (Forget & PieiTehumbert! T997[ ). Hence it was 
concluded that further climate simulations were required. 
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To investigate the possible climate of G1581d under a range 
of conditions, we have performed one-dimensional radiative- 
convective calculations. In Section |2] we discuss the model we 
used, while in Sectionl3]we present our results for varying atmo- 
spheric compositions, surface albedos, gravity and cloud cover- 
age. We also present some simple three-dimensional simulations 
that highlight the limitations of the one-dimensional globally av- 
eraged approach. In Section|4] we discuss the implications of our 
results and suggest directions for future research. 

2. Method 

Our radiative scheme is based on the coiTelated-A: model, with 
the absorption data calculated directly from high resolution spec- 
tra. We used the program kspectrunrj to produce line-by-Une 
spectra from the HITRAN 2008 database ( Rothmanet al.|2009| l. 
These were then converted to correlated-^ coefficients for use 
in the radiative transfer calculations. The CO2 collision-induced 
absorption (CIA) was parameterized using the results of lBaranov 
et aLJ ( |2004| l and [Gruszka & Borysow| ([19 98 ), while the sub- 
lorentzian profiles of iPerrin & Hartmann (11989} were used for 
the CO2 far line absorption. For further details and justification 
of the method used, refer to |Wordsworfli et ar] (|2010). Note that 
at high CO2 partial pressures (tens of bars), additional sources 
of CIA may appear that we have not been able to take into ac- 
count (Eymet et al. 2009 1. However, these would tend to increase 
the surface temperature in our calculations. As we are interested 
in making a conservative estimate of habitability, we can safely 
neglect them here. 

A matrix of coefficients was produced on a 6 x 9 temper- 
ature and log-pressure grid T - {100, 150, ... ,350) K, p = 

' http://code.google.eom/p/kspectrum/ 
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{10"\10"\...,10^) mbar. For mixed CO2-H2O atmospheres, 
we produced correlated-A: coefficients for H2O volume mixing 
ratios qn^o = {10^^ 10"^, . . . , 10"')- We used 32 spectral bands 
in the longwave and 36 in the shortwave. Sixteen points were 
used for the g-space integration, where g is the cumulated distri- 
bution function of the absorption data for each band. For the stel- 
lar spectra, we used the Virtual Planet Laboratory AD Leo and 
solar data ( Segura et al.|2005| l for M-class and G-class stars, re- 
spectively. The solar zenith angle 6^ was set to 60°, and the stellar 
flux was further divided by two to give a spatially and tempo- 
rally averaged value F,„ — 50/4 Vl - e^. As planetary albedo is 
a nonlinear function of zenith angle, the values obtained in ID 
can differ slightly from those produced by a three-dimensional 
model. However, a 60° zenith angle generally overestimates the 
globally averaged value, and hence predicts lower surface tem- 
peratures. We therefore neglected the error due to this difference 
here. Here Sq - 38L4W m"^ is the unaveraged stellar flux at the 
G1581d semi-major axis a, and e is the eccentricity. F„, is given 
along with a, e and the stellar luminosity L in Table L The factor 
involving e in the definition of F,„ comes from the fact that the 
flux averaged over one orbit is higher in an eccentric orbit than 
in an equivalent circular one. 

To account for the radiative effects of both clouds and 



Rayleigh scattering in our simulations, we used the Toon et al. 
(|1989 ) scheme. Rayleigh scattering was included by the method 



described in.Hansen & Travis 



( 1974 1, while the properties of the 



clouds were computed from Mie theory. Refractive indices for 
the Mie calculation were taken from ' Hansen| ( |2005l ) and | Warren] 
(|1984 ) for the CO2 and H2O clouds, respectively. In the simula- 
tions with CO2 clouds, a simple microphysical scheme was used 
to calculate the cloud opacity in each layer. Following [Forget 



et al. ( 2004 1, the number of condensable nuclei per kilogram of 



air Nc was treated as a free parameter. At each timestep, the mean 
cloud particle radius at each level was calculated from A^^. and the 
amount of condensed CO2. This was then used to calculate the 
mean particle sedimentation rates and cloud radiative properties. 
For water clouds, we used a simple approach. Fixed lay- 
ers of particles of radius 10 /^m and varying optical depth (de- 
fined at a reference wavelength L5 /vm) were placed in the low 
troposphere, which we defined arbitrarily as the level at which 
p — Q.5ps- In simulations with vertically varying amounts of 
water vapour, we used profiles of the form 



qH20 - qsatRH 



(1) 



where qHio is the water vapour mixing ratio, qsat the water 
vapour saturation ratio, and RH the relative humidity. The lat- 
ter was defined as in |Manabe & Wetherald| ( [T967| l 



RH = O.nip/psurf - 0.02)/0.98 



(2) 



until pipsurf < 0.02, after which we set qnio - 0. 

To model convection in the lower atmosphere we used ad- 
justment to keep the lapse rate F stable. In dry ideal atmospheres, 
r - -g/cp. As dense CO2 is a non-ideal gas, F can deviate 
from this value in the lower troposphere when the pressure is 
sufficiently high. We assess the importance of this deviation in 
the Appendix. The effects of moist convection on the lapse rate 
were not included. While moist lapse rates are shallower than 
dry ones, which results in lower surface temperatures, the quan- 
tity of unsaturated water vapour in cold atmospheres is low. The 
error due to this approximation was small in our calculations 
for surface temperatures below 273 K, and hence unimportant to 
our assessment of the minimum CO2 partial pressure necessary 
for habitability. In contrast, the inclusion of CO2 condensation 



Table 1. Planetary and stellar properties used in the one- 
dimensional simulations. When multiple values are shown, those 
in bold are the standard values used. 



Stellar luminosity 


^ [I^Sim] 


0.0135 


Stellar spectrum 




Sol, AD Leo 


Semi-major axis 


a[AU] 


0.22 


Eccentricity 


e 


0.38 


Averaged stellar flux 


F„[Wm-2] 


103.1 


Relative humidity 


RH 


0, variable, 1.0 


Surface gravity 


g[ms-'] 


10.0, 20.0, 30.0 


Surface albedo 


A, 


0.1,0.2,0.3 


CO2 partial pressure 


pco. [bars] 


0.2 - 50 


N2 partial pressure 


Pn, [bars] 


0, 1, 5, 20 


CO2 cloud cover 


Tl 


none, dynamic 


H2O cloud cover 


T2 


none, 1,2,3 



above the troposphere was extremely important. We took this 
effect into account using the algorithm described in Forget et al. 



(11998 ) and vapour pressure curves derived from Lide (2000). 

Thirty vertical levels in standard sigma coordinates were 
used, and all simulations were iterated in time until a steady 
state had been reached (this took up to 50 Earth years of sim- 
ulation time for the highest pressure runs). In all cases the en- 
ergy balance (incoming vs. outgoing radiation) of the model was 
checked at the end of each run. A list of all parameters used in 
the simulations is given in Table [T] 

3. Results 

We have studied the effects of a range of climate parameters on 
the mean surface temperature of G1581d. For clarity, this section 
is subdivided according to the effect studied. When not explic- 
itly stated otherwise in the text, the parameters used for each 
simulation are those in bold in Table 1 . 



3.1. Stellar spectrum 

To better understand the differences in climate caused by the fact 
that Gliese 58 1 is an M-class star, we first performed simulations 
comparing G-class (Sol) and M-class (AD Leo) stellar spectra. 
In both cases we normalized the total fluxes to the same value 
F„, (see Table [T]). Figure [T] shows the equilibrium mean surface 
temperatures obtained as a function of surface pressure for these 
two cases. 

As can be seen, the clear pure CO2 atmospheres under a G- 
class star collapse on the surface for pressures greater than about 
3 bar, but when the star is M-class, temperatures continue to in- 
crease, reaching the water melting point at just over 10 bar and 
a maximum value at around 30 bar. The essential reason for this 
difference is that Rayleigh scattering, which has an optical depth 
T« oc /l""*, has a much weaker effect on the red-shifted M-class 
stellar spectrum. 

In addition to reducing the effect of Rayleigh scattering, the 
red shift in the stellar spectrum leads to increased warming in 
the upper stratosphere. This effect is clear from the vertical tem- 
perature profiles, which are plotted in Figure [2] As can be seen, 
near-IR CO2 absorption leads to a temperature inversion near 
0.05 bar, which becomes even more pronounced if CO2 conden- 
sation is neglected in the calculation (Figurel2V)). 

To show the importance of these separate effects more 
clearly, in Figure l3] we have plotted the incident stellar spectra 
and the radiation reaching the ground for cases with and without 
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Fig. 1. Surface temperature vs. surface pressure for clear pure 
CO2 atmospheres under G-class (Sol) and M-class (AD Leo) 
stellar insolation. Dashed and dotted lines show CO2 and H2O 
phase transitions, respectively. In the G-class case, the line is 
discontinued because the atmosphere condenses on the surface 
for pressures greater than a few bar 



Rayleigh scattering included. As can be seen, Rayleigh scatter- 
ing dominates for wavelengths < 1 jum, while gaseous absorption 
dominates above 1 /i/m. 

As AD Leo has a mean temperature of around 3400 K, the 
planetary albedo values we obtain are lower than those for the 
blackbody 3700 K case that was used to define the habitable 
boundary in |Kasting et al.| ( [l993 i). This leads to a significant in- 
crease in the calculated surface temperature. Note that the tem- 
perature of G1581 is approximately 3200 K, so the planetary 
albedos we calculate using the AD Leo spectrum should lead 
to conservative habitability estimates. 

3.2. CO2 clouds 

If CO2 condensation occurs in the upper atmosphere of a planet, 
it should cause CO2 cloud formation. This effect is observed for 
example in the present-day Martian mesosphere (Montmessml 



et al. 20071. While there are still many unknowns associated 



with the microphysics of CO2 clouds, a basic estimate of their ef- 
fect can be gained through the scheme described in Sectionl2] In 
Figurel4k, the surface temperature vs. surface pressure is plotted 
for simulations with radiatively active CO2 clouds and varying 
values of the condensation parameter Nc- The infrared scattering 
effect described by Forget & Pierrehumbert (1997i is expected 



to be slightly less efficient when the star is M-class, because a 
greater portion of the incident stellar radiation is also scattered 
back to space. Nonetheless, the clouds still increase the surface 
warming by an amount that increases with Nc, up to a theoretical 
maximum of around 30 K at 40 bar for Nc = 10^ kg"^ FigureWp 
shows that at high pressures, the planetary albedo begins to in- 
crease significantly due to CO2 cloud coverage. This ultimately 
limits the maximum warming effect that the clouds can provide, 
as the downward scattering of thermal radiation from the planet's 
surface becomes less important than the direct reflection of stel- 
lar radiation to space. 

3.3. Surface gravity 




Fig. 2. Steady-state vertical temperature profiles for clear pure 
CO2 atmospheres of different surface pressures under M-class 
stellar insolation a) with and b) without CO2 condensation in- 
cluded. The dashed line shows the CO2 saturation pressure. 



Sotin et al. ( 2007 1 proposed a relationship between planetary 
mass and radius (r/rE) = aiM/ME'f, with a and/3 equal to LO 



and 0.274 for rocky planets and L262 and 0.275 for ocean plan- 
ets, respectively. Given this relation and the current uncertain- 
ties in G1581d's mass, we can expect the planet's surface gravity 
to be in the range 10 - 30 m s"^. In Figure p] we have plotted 
surface temperature vs. pressure for simulations with g - 10, 
20 and 30 m s"^. For a given atmospheric pressure, the CO2 
column amount pjg (and therefore the total mass of the atmo- 
sphere) decreases with g. Hence the primary effect of increasing 
g, predictably, is to cool the surface. 

Variations in g also change the adiabatic lapse rate, while 
leaving the CO2 saturation pressure unaffected. Hence for a 
given CO2 column amount the temperature profile varies with 
the gravity, which also influences the climate. These changes are 
important in determining the point at which maximum green- 
house warming occurs (see e.g. the g = 30 ms"^ case in Figure 
Bt), but they do not strongly affect the surface temperature for a 
given CO2 column amount at lower pressures. This can be seen 
in Figure |5]3, where the same three curves are plotted vs. pjg. 
The atmosphere begins to collapse before the surface tempera- 
ture reaches the water melting point only in the most conserva- 
tive ^ = 30 m s"^ case. 
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Fig. 3. Incident flux in the visible at the top of the atmosphere 
(black) and at the ground, with (blue) and without (red) Rayleigh 
scattering included. In this example the atmosphere consisted of 
40 bars pure CO2 (no clouds), while the stellar spectrum was 
G-class (Sol) for a) and M-class (AD Leo) for b). 
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Fig. 4. a) Surface temperature vs. surface pressure for pure CO2 
atmospheres under IVI-class insolation with CO2 cloud radiative 
effects included. Dashed and dotted lines show CO2 and H2O 
phase transitions, respectively, b) Planetary albedo vs. surface 
pressure for the same simulations as in a). 



3.4. Surface albedo 

Surface albedo has a dramatic efifect on the surface temperatures 
of planets with thin atmospheres, but its importance might be 
expected to be less in dense or cloudy atmospheres. We tested 
the effect of variations in surface albedo in our calculations for 
As - 0.1, 0.2 and 0.3. In Figures [6^ and b, we have plotted the 
surface temperature vs. pressure for clear-sky and A^c = 10^ kg"' 
CO2 cloud cases, respectively. 

As can be seen, the diff'erence in climate between A^ = 0.1 
and 0.3 is large in the clear-sky case. While the lower sur- 
face albedo leads to Ts > 273 K at 10 bar with a maximum 
Jj = 315 K at 40 bar, the higher surface albedo causes the at- 
mosphere to collapse on the surface for ps > 10 bar. However, 
in the cases with CO2 cloud cover (Figure l6b), the differences 
in surface temperature are generally smaller, and for all albedo 
values, surface temperatures above the water melting point are 
achieved before the surface pressure reaches 10 bar Hence if 
the atmosphere contains enough CO2 to increase surface tem- 
peratures above 273 K and clouds of some kind are present, the 
effect of surface albedo on the climate should not be critical. 



3.5. N2 partial pressure 

N2 is an major constituent of the atmospheres of several rocky 
bodies in the Solar System. We tested the effects it would have 
on the climate of G1581d by repeating the basic calculation of 
Section [TT] with added partial pressures of 1, 5 and 20 bars of 
N2. As there are no important N2 absorption bands in the visible 
or infrared, its primary effect on the climate is through a) ad- 
ditional Rayleigh scattering and b) changes in the specific heat 



capacity and mean molar mass of the atmosphere. We included 
both these effects in our calculations. N2 also has a warming ef- 
fect in the Earth's atmosphere due to the broadening of CO2 and 
H2O spectral lines (Goldbl att et al.|[20 09). However, this effect 
is relatively small (~5 K for Earth conditions) and we are inter- 
ested in conservative estimates of warming. Hence we neglected 
it here. 

As can be seen from Figure ITJ the effect of N2 on the sur- 
face temperature is relatively small. Only for p^^ = 20 bar is the 
planetary albedo significantly increased above the value for pure 
CO2 due to increased Rayleigh scattering. In this extreme case, 
the temperature fails to reach the water melting point, with the 
atmosphere collapsing after p^ = 20 bar However, if N2 warm- 
ing were included, this difference would presumably be less sig- 
nificant. 



3.6. H2O vapour and clouds 

Water plays an important and complex role in the radiative bud- 
get of Earth, through both direct vapour absorption and cloud 
absorption / scattering. Here we have studied its effects via the 
simple approach described in Section|2] Figure|8]shows runs per- 
formed with the relative humidity RH = at all levels, RH = 1 
at all levels, and RH defined as in (|2|. As can be seen, increas- 
ing water vapour has the predictable effect of increasing the to- 
tal warming, by an amount that increases with the surface tem- 
perature. This shifts the surface pressure at which we reach the 
water melting point down to 5 bar. Note that we neglect the ef- 
fects of moist convection on the tropospheric lapse rate, so the 
CO2+H2O results in Figure[8]at the highest pressures likely over- 
estimate the temperature somewhat. 
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Fig. 5. Surface temperature vs. surface pressure for clear pure 
CO2 atmospheres under M-class insolation with varying surface 
gravity. Surface temperatures are plotted vs. a) pressure and b) 
CO2 column density. Dashed and dotted lines in a) show CO2 
and H2O phase transitions, respectively. 
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Fig. 6. Surface temperature vs. surface pressure for pure CO2 
atmospheres a) without clouds and b) with CO2 clouds, Nc - 
10^ kg"', under M-class insolation with varying surface albedo. 
Dashed and dotted lines show CO2 and H2O phase transitions, 
respectively. 



Figure [9] shows simulations performed with fixed H2O 
cloud layers in the troposphere (see Section l2]l and varying va- 
por RH humidity as defined in (|2i. The planetary albedo in- 
creases with the cloud optical depth, but the amount by which 
it does so is lower at high pressures, because then atmo- 
spheric absorption/scattering above the clouds is more impor- 
tant. Counterintuitively, however, the clouds increase the surface 
temperature at low pressures and decrease it at high pressures. 
This effect can be explained by the fact that the atmosphere is 
relatively clear in some infrared spectral regions at low pres- 
sures. While an H2O cloud layer reflects starlight, increasing the 
planetary albedo, it will also absorb outgoing infrared radiation 
from the surface, eff'ectively acting as a 'thermal blanket' and in- 
creasing the surface temperature. At high pressures, however, the 
lower atmosphere absorbs strongly in almost all infrared bands 
due to increased water vapour and the effects of CO2 CIA. Hence 
the clouds have a minimal effect on the outgoing infrared radia- 
tion, but still increase the planetary albedo, which leads to a net 
cooling. 

For higher values of H2O cloud optical depth, we found 
that at intermediate pressures, the atmosphere collapsed on the 
surface (see also Figure 10 1. This suggests that sufficient H2O 
cloud cover at low altitudes may limit the habitability of G1581d. 
However, we used 100 % fractional cloud coverage in these sim- 
ulations, while on Earth, global low-level H2O cloud coverage is 
generally less than 50 %, with reported optical depths of 4.7 at 
0.6 //m (Kitzmann et al. ( 2010|l, deriv ed from results of'Rossow' 
^chiffer (1999) and [Warren et al.| ( |2007 )). Hence our results 
for large optical depths may overestimate the cooling effect of 
the clouds at high pressures. For a more comprehensive study 



of the effects of H2O clouds than we have given here, three- 
dimensional climate simulations including cloud microphysics 
are necessary. 

3.7. Most conservative scenario 

What happens if we combine the most conservative estimates for 
surface gravity, surface albedo, N2 abundance and CO2 cloud 



cover? Figure 10 shows a simulation that was performed with 

g = 30 m s"^ 
erage (Nc = 



As = 0.3, pn^ - 20 bar, minimal CO2 cloud cov- 
10^ kg"') and water vapour present according to ^. 
The H2O cloud opacity was varied between t-1 and 3. As can 
be seen, in all three cases the atmosphere collapses at CO2 pres- 
sures between 3 and 20 bar. Nonetheless, it can be seen that at 
higher pressures, the atmosphere is once again stable, with sur- 
face temperatures well above the water melting point at pcoi = 
40 bar for these values of the H2O cloud opacity. Given the initial 
assumption of CO2 partial pressure of tens of bars, therefore, the 
outlook for the habitability of G1581d would appear rather good. 

3.8. Tliree-dimensional equilibrium temperatures: airless 
planet case 

While one-dimensional simulations of the type described here 
are useful because of their simplicity, they make the inevitable 
assumption of globally averaged stellar insolation. As discussed 
in S elsis et al. (200 7), G1581d is likely to be in a tidally resonant 
orbit, which means that there will be dramatic differences in in- 
solation across the surface of the planet. To get an idea of the 
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Fig. 7. Surface temperature vs. surface pressure for clear mixed 
CO2-N2 atmospheres under M-class insolation. Dashed and dot- 
ted lines show CO2 and H2O phase transitions, respectively, b) 
Planetary albedo vs. surface pressure for the same simulations as 
in a). In these figures, the jc-axis shows the CO2 partial pressure. 
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Fig. 8. Surface temperature vs. surface pressure for clear mixed 
CO2-H2O atmospheres under M-class insolation. Dashed and 
dotted lines show CO2 and H2O phase transitions, respectively. 



effect this could have on the climate, we performed some simple 
equilibrium temperature calculations in three dimensions. 

In Figure [TT] we have plotted the temporal evolution of sur- 
face temperature on the 0° and 180° longitude sides of an airless 
planet with the same orbit as G1581d, zero obliquity, tidal reso- 
nances of 1:1 (left) and 2:1 (right), surface albedo A, = 0.2 and 



Fig. 9. a) Surface temperature vs. surface pressure for mixed 
CO2-H2O atmospheres under M-class insolation with no CO2 
clouds and fixed H2O cloud coverage. Dashed and dotted lines 
show CO2 and H2O phase transitions, respectively, b) Planetary 
albedo vs. surface pressure for the same simulations as in a). In 
all cases optical depth is defined at the reference wavelength (1.5 
//m). 

initial uniform surface temperature 250 K. Note that the 'dark 
side' temperature in the plot on the left (solid line) is still de- 
creasing toward K in the simulation due to the finite thermal 
inertia that was given to the surface. 

The dotted lines in the figure correspond to the equilibrium 
sub-stellar temperature at periastron 
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As can be seen, at one location the local surface temperature can 
vary by as much as 300 K over the course of one orbit. This 
variability reduces for planets with atmospheres as a function of 
the surface pressure - on Venus, for example, where the diurnal 
period is 171 Earth days and the pressure is 93 bars, the only sig- 
nificant surface temperature variation occurs due to topographic 
effects ( Mueller et al.|2008| . However, the decrease in variability 
of r, with pressure will be dependent in general on the specific 
atmospheric and orbital parameters of each planet. Hence fur- 
ther modeling is required before global temperature variations 
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Fig. 10. Surface temperature vs. surface pressure for the 'most 
conservative case' described in the text. Dashed and dotted 
lines show CO2 and H2O phase transitions, respectively, while 
crosses, circles and stars indicate surface temperatures for atmo- 
spheres with H2O clouds of opacity t = 1, 2 and 3, respectively. 
In this figure, the ;c-axis shows the CO2 partial pressure. 



can be estimated for the hypothetical G1581d atmospheres we 
have studied here. 

The second point of interest raised by this simple 3D model 
is that locally, the temperature on the planet rises well above the 
water melting point during the 'warm season' (substellar point at 
periastron), despite the fact that there is no atmospheric warm- 
ing. This raises an interesting question: could a mainly frozen 
planet in this orbit with a thin (p, > 6.1 mbar) atmosphere con- 
tain a zone where liquid water forms seasonally? Two extreme 
cases can be envisaged, depending on the total amount of wa- 
ter present. If the fraction of water on the planet is small, the 
regions of low solar insolation are likely to act as cold traps in 
which water ice forms, leaving the zone in which liquid water is 
possible dry. If the planet's surface is covered by H2O, however, 
there will be no limit to the amount of local evaporation that can 
occur. In this case, the local climate could potentially be quite 
interesting, although the high albedo of water ice would reduce 



the surface temperatures seen in Figure 1 1 



4. Discussion 

Through one-dimensional radiative-convective climate model- 
ing, we have found that as little as 5 bars of CO2 (~2.5 bar equiv- 
alent column amount under Earth gravity) may be sufficient to 
maintain a global mean temperature above the melting point of 
water. For CO2 partial pressures greater than 40 bar (~13 bar 
equivalent column amount under Earth gravity), we found that 
the mean surface temperature of Gliese 58 Id is high enough to 
allow surface liquid water even when the most conservative val- 
ues of gravity, surface albedo and N2 partial pressure are chosen, 
although strong dependence on the total H2O cloud opacity was 
observed. 

We tested the dependence of surface temperature on a wide 
range of atmospheric and geophysical parameters. Because the 
M-class spectrum of Gliese 581 is red-shifted with respect to 
that of the Sun, heating due to near-IR absorption in the upper 
atmosphere increases in our simulations, while Rayleigh scatter- 
ing is much less important. This is in agreement with the general 
habitability study of Kasting et al. ( 1993 1, but we find the dif- 
ference to be even stronger than was reported there, due to the 
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Fig. 11. Snapshot surface temperature on the points of a hypo- 
thetical tidally resonant Gliese 58 Id with minimal atmosphere 
and surface albedo A - 0.2 that are sub- (dashed line) and anti- 
stellar (solid line) at periastron. Upper and lower dotted lines in- 
dicate the equilibrium sub-stellar temperature at periastron and 
apastron, respectively. The initial surface temperature in the sim- 
ulations was 250 K at time - days. 



low effective temperature of G1581. We also estimate an mod- 
est warming effect due to CO2 clouds (up to 30 K at pco2 - 50 
bar), dependent on details of the cloud microphysics, although 
the efficiency of this process is reduced compared to the case for 
G-class stars. 

We find that increasing gravity primarily decreases the CO2 
column amount at which the maximum surface temperature 
occurs. Unsurprisingly, increased surface albedo decreases the 
mean surface temperature, although if there is cloud coverage, 
this effect is relatively small for the high CO2 pressures at which 
the surface starts to become habitable. Increasing the partial 
pressure of N2 has little effect below partial pressures pt^^ of tens 
of bars. Finally, water vapour increases the surface temperature 
by a moderate amount below Tj = 273 K, after which its role 
becomes increasingly important. We briefly assessed the impact 
of H2O clouds, and found that they warmed (cooled) the surface 
at low (high) pressures, due to the interplay between their effects 
on incoming stellar and outgoing thermal radiation. 

As discussed by [Selsis et al] ( |2007| l, N2-CO2-H2O atmo- 
spheres likely only represent a fraction of the possible range 
for Earth-like / super-Earth extrasolar planets. While the pres- 
ence of radiatively inactive gases such as argon would be of lit- 
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tie importance to the calculations discussed here, small amounts 
of other greenhouse gases such as methane would result in in- 
creased absorption, if such gases could remain chemically stable 
in the atmosphere. Unless the near-IR absorption of these gases 
were large enough to cause a strong temperature inversion in the 
upper atmosphere, this would further increase the surface tem- 
perature of the planet. The presence of other clouds or haze lay- 
ers could alter the radiative balance in less predictable ways, de- 
pending on their scattering / absorption properties and height in 
the atmosphere. Unfortunately, effects of this kind are difficult to 
constrain further until direct atmospheric observations become 
available. 

The results discussed here clearly have interesting implica- 
tions, as they show that G158 Id could be the first discovered hab- 
itable exoplanet. They are also testable by future observations - 
dense CO2 atmospheres have recognizable spectral signatures 
that could be detected, for example, using the proposed Darwin 
or TPF missions. However, there is an important limitation to 
this study: all simulations performed were one-dimensional. As 
we noted in Section |3.8| calculations of the globally averaged 
surface temperature neglect variations due to local changes in 
stellar insolation (as well as topography and other effects). This 
has serious implications for dense, relatively cold CO2 atmo- 
spheres, as if the dark side temperature of the planet is too low, 
CO2 will condense on the surface. To test whether a dense atmo- 
sphere on G1518d would be stable against the collapse of CO2 in 
this way, we plan to repeat the calculations reported here in the 
future using a three-dimensional climate model. 
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5. Appendix 

To assess the importance of the non-ideal behaviour of CO2 
gas at high pressures in our calculations, we used the modified 
Redlich-Kwong equation of state proposed by Li et al. (2006i. 
The Redlich-Kwong equation takes the form 



RT 



v-b v{v + b)^lf 



(5) 



where p, T, v and R are pressure, temperature, specific volume 
and the ideal gas constant, respectively. Various definitions of 
a and b are possible; in Li et al. (2006 1 a - ay(T) + a2{T)p, 
with a I and 02 polynomial functions of temperature, while b = 
0.0E664RTc/pc, with T^ and pc the critical point temperature and 
pressure for CO2. Equation (51) can be written as a cubic equation 
in terms of v and hence is solvable analytically. The inverse com- 
pressibility factor/? = 2"' - RT/pv calculated from ([Sjl agrees 
approximately with the values stated in Kasting (1991) for the 
triple point and saturated CO2 at 0°C. This is unsurprising, as 
the equation used there was of the same basic form. 

To get a first-order estimate of the effect of non-ideality on 
the tropospheric lapse rate (and hence on the climate), we ap- 
proximated (|5]) to a virial form as 



pv/RT ^l+B/v 

with 

B(T) = fe ■ 



RTV2- 



(6) 



(7) 



Given the definition of specific heat capacity in Kasting 
(|T99T) 



cpip, T) = c,(0, T)-T 



X 






pdp 



it is clear that in the virial approximation, to leading order 
Cpip, T) = Cp(Q, T) - TpB"{T) 



(8) 



(9) 



Various formulations for Cp(0, T) were considered, including 
the quadratic f unction given in [K asting (1991) and the power 
law function of [Lebonnois et aL[ ( ,2010j . We found that the lat- 
ter expression, which is defined as Cp{0,T) = 1000(^/4 60)'"^ 
had the greatest accuracy vs. reference data (Lide"2000). Using 
this and calculating B"(T) numerically by finite differences, we 
found for T = 250 K and p = 0, 50 bar that Cp = 795.9 and 
825.0 J K"' kg"^ respectively. Climate calculations performed 
with fixed specific heat capacities of 750 and 850 J K"' kg"' 
showed maximum temperature differences of only 2.2 K over 
the entire range 0.2 - 50 bar, with the latter value of Cp giving 
colder results. Hence we decided to ignore variations in Cp with 
altitude in our calculations, and used the upper limit Cp = 850 
J K ' kg ' everywhere in the pure CO2 calculations. For mixed 
CO2-N2 atmospheres, a weighted mean value for Cp was used, 
with the heat capacity of N2 taken to be a constant 1003 J K~' 
kg-'. 



